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ABSTRACT 


The  proh.'.em  of  supersonic  and  hypersonic  flow  past  blunt  bodies  at 
small  angles  of  attack  is  considered.  Two-dimensional  asymmetric  as 
well  as  three-dimensional  flow  is  analyzed. 

The  method  of  analysis  is  an  inverse  one^  that  is^  the  shock-wave 
shape  and  free  stream  conditions  are  known^  and  the  corresponding  body 
shape  and  flow  field  are  to  be  determined.  A  semi-analytical^  semi- 
numerical  solution  to  the  full  inviscid  equations  of  motion  is  obtained 
by  expanding  the  stream  function  and  density  in  double  McLaurin  aeries 
in  powers  of  distance  from  the  shock-wave  axis  of  symmetry  and  powers  of 
the  shock-wave  angle  of  attack.  These  expansions  bring  about  a  separa¬ 
tion  of  variables  that  reduces  the  problem  to  numerical  integration  of 
ordinary  differential  equations  for  the  power-series  coefficients.  The 
integration  can  be  accomplished  only  by^truncating  the  series.  This 
truncation  may  be  performed  after  keeping  any  desired  number  of  terms 
in  the  series  about  the  shock  symmetry -axis .  Only  first-order  terms 
in  shock-wave  angle  of  attack  are  retained. 

Solutions  at  zero  angle  of  attack  are  obtained  as  a  special  case 
of  the  general  problem.  Results  at  zero  angle  are  obtained  for  shock 
waves  that  are  portions  of  circles,  parabolas,  spheres,  and  paraboloids 
of  revolution  at  a  free  stream  Mach  number,  M  ,  of  infinity  and  ratio 
of  specific  heats,  j  ,  of  l.U.  The  first  three  truncations  are  carried 
out  for  the  parabolic  and  spherical  shock,  whereas  the  first  four 
truncations  are  calculated  for  the  circular  and  paraboloidal  shock.  In 
all  cases  the  stagnation-point  density  is  obtained  to  five-significant- 
figure  accuracy  in  two  truncations.  Standoff  distance  is  obtained  to 
four -place  accuracy  in  three  truncations.  Convergence  away  from  the 
axis  of  symmetry  is  not  so  rapid.  Body  shapes  out  to  the  sonic  point 
converge  in  three  truncations  for  the  two-dimensional  cases.  However, 
axisjnrmetric  body  shapes  and  both  t-wo -dimensional  and  axisymmetric  sonic 
line  positions  and  surface  pressure  distributions  out  to  the  sonic  point 
are  not  converged  in  four  truncations,  but  are  very  close  to  being  so. 


ii  - 


Comparison  with  results  obtained  using  Van  Dyke's  and  Garabedian's 
numerical  solutions  indicates  that  the  method  under  consideration  is 
more  accui'ate  than  the  Van  Dyke  method. 

The  first  three  truncations  are  calculated  for  parabolic  and  parab¬ 
oloidal  shock  waves  at  small  angle  of  attack  and  infinite  free  stream 
Mach  number.  The  specif ic -heat  ratio,  y  ,  takes  on  the  vtilues  l.U, 

1.2,  1.1,  and  1,05.  For  all  cases,  the  streamline  that  wets  the  body 
passes  through  the  shock  wave  slightly  above  the  point  where  the  shock 
wave  is  normal,  and  consequently  does  not  possess  maximutri  entropy. 
Although  these  results  provide  a  counter  example  to  the  conjecture  that 
any  isolated  convex  body  in  a  supersonic  stream  is  wetted  by  the  stream¬ 
line  of  maximum  entropy,  the  stagnation  and  maximum -entropy  streamlines 
are  very  close  to  each  other,  thus  suggesting  that  the  conjecture  may 
be  a  good  practical  approximation.  Indications  are  that,  in  the 
Newtonian,  limit  (M=oo^  7=1),  the  spnjecture  is  true  for  the 
shock  shapes  considered- 

A  symmetrical  shock  wave  at  angle  of  atta.ck  is  produced  by  a  body 
that  is  asymmetric  to  the  shock  axis  of  symmetry.  However,  a  conic 
section  may  be  fitted  to  the  converged  body  shape  out  to  the  sonic  point. 
This  is  done  for  the  third-truncation  body  that  produces  a  parabolic 
shock  at  angle  of  attack  of  10°.  The  body  is  closely  approximated  by 
a  prolate  ellipse  at  angle  of  attack  of  1U.2°.  Thus,  for  the  two- 
dimensional  case,  the  shock  wave  is  seen  to  rotate  more  slowly  than 
the  body  as  angle  of  attack  is  imposed.  The  same  result  is  true  for 
the  third-truncation  body  in  the  three-dimensional  case. 
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I.  INTRODUCTION 


For  more  than  a  decade ,  theoreticians  interested  in  high-speed 

gasdynaraics  have  been  seeking  a  solution  to  the  problem  of  supersonic 

and  hypersonic  flow  past  a  blunt  body  that  combines  the  advantages  of 

minimum  computational  difficulty  with  maximum  accuracy  of  results. 

Physically,  the  problem  under  consideration  is  the  following.  A  blunt- 

nosed  configuration  is  traveling  through  a  uniform  gas  at  a  flight  Mach 

number  greater  than  unity.  A  detached  shock  wave  forms  ahead  of  the 

body  at  a  position  depending  on  the  body  shape,  flight  direction  and 

Mach  number,  and  certain  physical  properties  of  the  gas.  The  gas  ahead 

of  the  shock  wave  is  unaffected  by  the  presence  of  the  body.  It  is 

desired  to  determine  the  flow  properties  (e.g. ,  velocity,  pressure, 

density)  in  the  region  between  the  shock  and  the  body.  In  general, 

this  region  may  be  divided  into  two  subregions,  depending  on  the  magni- 

tude  of  the  local  flow  Mach  number  with  respect  to  unity.  The  locus  of 

points  at  which  the  local  flow  Mach  number  is  unity  is  called  the  sonic 

line.  If  one  considers  steady  flow  relative  to  the  body  and  denotes 

■* 

by  M  and  M^  the  free  stream  and  local  Mach  numbers  respectively, 
then  an  example  of  the  above  situation  is  here  depicted. 
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Mathematically^  the  problem  considered  is  governed  by  highly  non¬ 
linear  differential  equations  in  a  domain  possessing  a  free  boundary 
along  vhich  initial  values  may  be  determined  if  the  boundary  shape  and 
free  stream  conditions  are  known„  The  governing  equations  are  of 
elliptic  type  in  the  subsonic  region  (  <  1  )  and  of  hyperbolic  type 

in  the  supersonic  region  (  >  1  ) .  The  initial-value  problem  for 

hyperbolic  equations  is  well  posed  and  yields  to  solution  using  the 
well-known  theory  of  characteristics.  The  initial-value  problem  for 
elliptic  equations  is,  however,  poorly  posed,  and  much  of  the  effort 
expended  on  the  blunt-body  problem  has  been  confined  to  its  solution  in 
the  subsonic  region.  Attempts  have  been  made  to  solve  the  problem  by 
both  inverse  and  direct  methods.  In  the  inverse  m.ethod,  the  shape  of 
the  detached  shock  and  the  free  stream  conditions  are  prescribed,  and 
the  corresponding  body  shape  and  flow  field  are  to  be  determined.  In 
the  direct  method,  the  body  shape  and  free  stream  conditions  are  pre- 
scribed,  and  the  corresponding  shock  wave  shape  and  flow  field  are  to 
be  determined.  References  [1]  and  [2]  contain  detailed  reviews  of 
and  extensive  references  to  previous  work  for  plane -symmetric  and  axi- 
symmetric  shapes  at  zero  angle  of  attack,  i.e.,  bodies  whose  flight 
direction  is  parallel  to  the  body*  axis  of  symmetry. 

The  methods  now  in  popular  use  for  solving  the  zero -angle  problem 
in  the  entire  subsonic  region  to  a  high  degree  of  accuracy  are  numerical 
in  nature  and  require  the  use  of  high-speed  electronic  computing  machine 
Two  such  methods  are  Van  Dyke's  [3]  approach  to  the  inverse  problem  and 
Belotserkovskii ' s  [4]  solution  to  the  direct  problem.  Van  Dyke's 
analysis  essentially  reduces  the  solution  of  the  conservation  equations 
for  steady,  inviscid  flow  of  a  perfect,  non-heat-conducting  gas  to 
numerical  integration  of  two  coupled  partial  differential  equations  for 
the  density  and  stream  function.  Belotserkovskii,  on  the  other  hand,  by 
assuming  a  polynomial  variation  between  the  shock  and  body  for  certain 
flow  quantities,  reduces  the  problem  to  integration  of  ordinary  differ¬ 
ential  equations.  Practical  experience  has  shown  that  both  of  these 
methods  require  a  large  amount  of  skill  on  the  part  of  the  numerical 
analyst  in  order  to  obtain  satisfactory  solutions  for  other  than  simple 


"body  shapes  such  as  spheres  or  circular  cylinders.  Hence,  even  though 
these  methods  yield  fairly  high  accuracy,  they  possess  a  moderate  amount 
of  computational  difficulty. 

In  contrast  to  the  zero-angle  problem,  comparatively  little  effort 
has  been  expended  on  the  problem  of  the  blunt  body  at  angle  of  attack 
in  either  the  two-dimensional  or  three-dimensional  cases.  Both  Mangier 
[5]  and  Vaglio-Laurin  and  Ferri  [6]  suggest  computational  methods  for 
attacking  such  problems,  but  neither  presents  results  using  the  analyses 
proposed.  They  both  obtain  solutions  to  the  zero -angle  inverse  problem 
by  introducing  the  stream  function  as  an  independent  variable  through  a 
von  Mises  transformation.  The  resulting  partial  differential  equations 
are  then  numerically  integrated  from  the  shock  to  the  body  using  a 
marching  technique.  Mangier  suggests  an  extension  of  this  marching 
technique  to  the  three-dimensional  angle -of -attack  problem  formulated 
in  terms  of  a  pair  of  stream  functions.  THe  approach  of  Vaglio-Laurin 
and  Ferri  to  the  problem  of  the  blunt  body  at  incidence  involves  a 
first-order  perturbation  of  the  flow  variables  about  their  zero-angle 
values. 

Vaglio-Laurin  [7]  has  developed  another  approach  to  the  blunt -body 
problem  at  zero  or  finite  incidence  which  involves  the  successive 
refinement  of  an  approximate  solution  by  application  of  the  PLK  method. 
Results  for  configurations  at  both  zero  and  finite  angle  of  attack 
using  the  analysis  are  presented  in  ref.[T]. 

The  success,  however,  of  Mangler's  and  Vaglio-Laurin ’ s  methods  for 
the  solution  of  the  angle -of -attack  problem  depends  on  an  assumption 
regarding  a  basic  property  of  asymmetric  supersonic  flows.  The  assump¬ 
tion  is  that  the  body  is  covered  by  the  streamline  that  crosses  the 
shock  at  right  angles  and  thus  possesses  maximum  entropy.  Mangier  has 
conjectured  that  such  a  condition  might  be  a  feature  of  any  flow, 
whereas  Vaglio-Ijaurin  [7,  P-19]  purportedly  has  proven  this  to  be  true 
for  two-dimensional  asymmetric  configurations.  However,  the  proof  has 
been  found  to  be  faulty,  and  Vaglio-Laurin  is  in  accord  with  this  con¬ 
clusion  (private  communication). 
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The  analysis  that  follows  yields  solutions  for  two-and  three- 
dimensional  supersonic  and  hypersonic  flow  past  hlunt  "bodies  at  small 
angles  of  incidence.  The  method  is  an  inverse  one,  that  is,  the  shock 
wave  shape,  free  stream  conditions,  and  angle  of  incidence  of  the  shock 
wave  are  known,  and  the  corresponding  body  and  flow  field  are  to  be 
obtained.  The  equations  of  motion  are  solved  using  a  semi-analytical, 
semi -numerical  approach  which  reduces  the  problem  to  numerical  Integra 
tion  of  ordinary  differential  equations.  The  solution  may  be  carried 
to  any  desired  degree  of  accuracy  and  requires  no  extensive  numerical 
analytical  skill.  Moreover,  no  assumption  need  be  :.iade  regarding  the 
behavior  of  the  maximum-entropy  streamline.  Its  position  in  the  flow 
is  obtained  as  part  of  the  solution.  Results  at  zero  incidence  arc 
obtained  as  a  special  case  of  the  general  problem.  Although  the  exam¬ 
ples  chosen  are  for  specific  shock  shapes  at  infinite  Mach  number  and, 
in  the  majority  of  cases,  a  specif ic -heat  ratio  of  l.k,  the  analysis 
is  performed  for  the  general  case.  The  niethod  of  extension  to  higher 
angles  is  obvious  and  introduces  no  new  difficulties. 
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II.  ANALYSIS 


A.  Coordinate  System 

The  coordinate  system  adopted  is  that  developed  by  Van  Dyke  [3]. 

The  detached  shock  wave  is  described  by  a  conic  section  so  that  it  is 
a  portion  of  a  hyperboloid  of  revolution,  paraboloid,  prolate  ellipsoid, 
sphere,  or  oblate  ellipsoid.  In  cylindrical  polar  coordinates  originat¬ 
ing  from  its  vertex  (Pig.  l) ,  any  such  shock  has  the  equation 

r^  =  2R  X  -  Bx^  (l) 

r  '  ' 


where  R  is  the  nose  radius  of  the  shock  and  B  is  a  parameter  char- 
s 

acterizing  the  eccentricity  of  the  conic  section.  B  is  negative  for 
hyperboloids,  zero  for  paraboloids,  positive  for  ellipsoids,  and  unity 
for  a  sphere.  >» 

An  orthogonal  coordinate  system  (  |,ti  )  that  contains  the  shock 
wave  as  a  coordinate  surface  is  defined  by 


x 


R 


5 


/(I  -  B|^)(l  -  B  +  Bti^) 

J 


(2a) 


^  (2b) 

s 

Special  cases  are 

^  =  ^(l  +  1^-  rj^)  for  the  parabola  (B  =  O)  (2c) 

s 


^  =  1  -  ri  \/l  -  for  the  circle  (B  =  l)  (2d) 

n 

s 


The  shock  wave  is  described  by  q  =  1  . 

In  this  coordinate  system,  with  the  azimuthal  angle  <p  as  the 
third  orthogonal  coordinate,  the  differential  line  element  ds  is  given 
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(3) 


by 


,  d.  »  2  a. 2  ,  ?  ,  2  ,  ^  ,  2 

dc  ==  d|  +  ai]  +  h^  d</j 


/ 


1  -  Bl 


&i  + 


.^2  ^ 
4  +  n 

c  +  Bn'" 


dT|  +  v|  Ti‘'  d(p 


whei-e  C  s  1  -  B  .  Here  v  is  a  geometric  parameter,  being  0  for 
plane  and  1  for  axisymmctric  or  three-dimeneioiial  flow. 


B ,  Differential  Equations,  Initial  Conditions,  and  Method  of  Solution 

Denote  by  q,  .o,  and  p  the  velocity,  density,  and  pressure  nori- 

diitienaionnlized  with  respect  to  the  free  stream  speed  ,  the  free 

stream  density  p  ,  and  twice  the  free  stx’eam  dynamic  pressure 
?  “ 

ly  ■  Then,  the  equations  expressing  conservation  of 
mass,  momentum,  and  entropy  for  steady  inviscld  flow  of  a  perfect,  non¬ 
heat  -conducting  gas  having  constant  specific  lieats  are 


V.(o?)  =  0 

{ha) 

(iq  .Vq  -1-  Vp  =  0 
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the  adiabatic  oxiionent. 
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Hero  u,v  ,  and  w  are  the  components  of  q  in  the  I ,  >  and  <p 

dii'ccLlcnn,  and  subscripts  i,  q  ,  and  <p  denote  partial  differentla- 
iioi]  in  tile  usual  sense. 

eince  the  method  of  solution  varies  silghtly  between  the  twc-and 
;  iin  e -d  imcns  iona  1  cases,  departiire  will  now  be  taker;  from  •:  c^-rcur  r^nt 
;  rc’ at  inert . 


1.  Plane  Flow 


Equations  (5a)  -  (5e)  apply  with  first 

(continuity)  equation  is  satisfied  by  introducing  a  stream  function 
defined  by 


2  2 
eg  +  n 


%  =  v/  ^  PU 

‘  V  c  +  Bq 


■^g  = 


/ 


1  -  Bg^ 


(6) 


then  the  last  (entropy)  equation  simply  states  that 


p  =  p^f(t) 


(7) 


where  the  functiojial  form  of  f  is  determined  by  the  boundary  conditions 
at  the  shock.  Using  equation  (?)  to  eliminate  the  pressure  from  the 
equations  of  motion  yields 
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Thus  far,  the  attitude  of  the  body  with  respect  to  the  free  stream, 
or,  since  we  are  considering  the  Inverse  problem,  the  angle  of  the 
uniform  free  stream  with  the  shock  axis  of  symmetry,  has  not  entered 
Into  the  mathematical  formulation.  The  angle  of  attack  of  the  shock 
will,  however,  appear  in  t.ne  boundary  conditions  and  in  the  explicit 
expression  for  f('^)  .  These  conditions  will  be  derived  by  employing 
McLaurin  series  expansions  about  zero  angle  of  attack  and  retaining 
first-order  terms  in  shock  angle  of  attack,  e  .  Setting  6  equal  to 
zero  recovers  the  zero-angle  conditions. 


Moreover,  in  the  angle-of -attack  problem,  the  position  in  the  flow 
ahead  of  the  shock  wave  of  the  streamline  that  wets  the  body  is  not 
known  a  priori ,  but  must  be  obtained  as  part  of  the  solution.  A  funda¬ 
mental  question  in  the  analysis  of  hypersonic  flows  past  inclined  blunt 
bodies  concerns  the  behavior  of  this  streamline.  It  has  been  conjectured 
[5,7]  that  any  isolated  convex  body  in  a  uniform  supersonic  stream  is 
covered  by  the  fluid  of  maximum  entropy  In  the  now,  This  means  that 
the  stagnation  streamline  must  have  crossed  the  bow  shock  where  it  is 
a  normal  shock  wave.  To  determine  the  position  of  the  stagnation  stream¬ 
line  in  the  free  stream  and  thus  possibly  disprove  the  conjecture,  the 
I’ormuiation  of  tlio  angle-of -attack  problem  requires  the  Introduction  of 
an  additional  parameter,  6  (Fig.  2).  The  quantity  cb  is  a  measure 
ol'  the  distance  ol'  the  body  streamline  in  the  free  stream  from  a  given 
reference  line.  This  reference  line  is  at  angle  g  to  the  x-axls  and, 

X  1  r 

in  general,  is  taken  to  pass  through  the  point  ~  “  5  ^  ~  ~  ^  ' 
llowc'vc'r,  in  the  singular  case  of  the  parabola  (  B  =  0  ),  the  refcrcrife 
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Values 
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u  I  j'cc  t  icn . 


en  through  the  focus  of  the  parabola  (  ^  ~  ^  ^  W~  ~  ^  ^ ' 

of  u,  V,  p,  |j  just  behind  the  shock  wave  (  t]  =  1  )  are 
the  oblique  shock  relations  (c.g.,  ref.  [3])  in  terms  of  the 
the  tangent  to  the  bow  wave  makes  with  the  free  stream 
Expressing  this  angle  in  terms  of  the  slope  of  the  bow  wave 
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using  eqr. .  (n),  and  retaining  first  order  terms  in  0 
ai  conditions 
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at  n  =  1 


and  gives  for  the  function  f('>Jr) 
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whore  M  is  the  free  c-troam  Mach  number.  For  the  .singu.lar  case  of  a 
])ai-aboiic  shock  wave  (  B  -  0,  C  =  .1  ),  the  ex])re3sions  for  p,  i|f,  ij/  nt 
1  =  1  j  .'U’.d  lor  f  are 
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The  boundary  condition  at  the  tody  is 

\Jr  =  0  at  Ti  =  ^  (11) 

The  method  of  solution  essentially  consists  of  separation  of 
variables  in  equations  (8)  -  (ll)  by  expanding  the  dependent  variables 
in  double  McLaurin  series  about  the  shock  axis  of  symmetry,  i  =  0  , 
and  about  e  =  0  .  Only  first  order  terms  in  e  are  retained.  This 
leads  to  ordinary  differential  equations  with  corresponding  boundaiy 
conditions  for  unkiiown  functions  of  q  which  can  bo  soJ.ved  only  by 
truncating  the  series.  The  resulting  equations  are  readily  integrated 
by  standai'd  numerical  methods.  The  above  procedin’c  will  now  be  carried 
out . 

Set 


+  €[^'^o(h)  +  f,j(q)i"+  o(^^^)l  +  o(€^)  (12:.) 
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(  1  ^b ) 


T'he  first  subscript  on  the  I'unctions  of  is  0  for  the  zero  angle  of 

attack  nortion  of  the  solution  and  1  for  the  angle  of  attack  portion, 
and  the  second  subscript  indicates  the  position  of  the  term  in  its 
respective  series.  The  forms  of  these  expansions  are  readily  ascertained 
by  consideration  of  the  initial  conditions. 

Expansions  (l2a)  and  (l2b)  are  substituted  into  the  differential 
equations  (8)  and  boundary  conditions  (9)  -  (ll)*  Equating  coefficients 
of  like-  powers  of  5  and  t  results  in  successive  problems  for 

SqO'  Soi  ^  "'oo^  ^00>  ^01>  ^0i>-S02  ’  zero-order  problem 

for  f__  and  involves  two  equations  which  also  conta.ln  g--,  ,  the 

00  00  Oj. 

first-order  problem  for  ^qq>  ^oi’  ir'^olves  four  equations 

which  also  contain  g^^  ,  and  so  forth.  The  same  ivittern  occurs  in  the 
successive  equations  for  the  angle-of -attack  functions.  This  "backward 
influence"  of  n^^-order  problem  is  due  to  the  ellipticity 

of  the  differential  equations  in  the  subsonic  region^  and  is  a  mathc- 
m<atlcal  manifestation  of  the  physical  fact  tiiat  disturbances  are  propa¬ 
gated  in  ail  directions  in  a  subsonic  flow.  This  mathenxitica]  behavior 
may  be  contrasted,  for  oxami'le,  with,  tliat  of  the  equations  determirung 
the  coefficients  of  the  Biasius-series  solution  to  th,c  parabolic 
cqu.atiuns  governing  the  laminar  boundary  layer  on  n.  blunt  body,  .lub- 
stitution  of  the  Biasius  scries  into  the  governing  dil’ferential  equation 
and  equ.ating  coefficients  (jf  like  powers  result.s  in  (>quations  for  the 
sericr,  coefficients  wliicli  involve  only  one  coefficient  or  depend  only 
on  Jower  order  coefficients,  and  hence  may  be  solved  success iveJ.y . 

t  h 

We  return  to  the  problem  under  consideration.  To  solve  the  n  '- 
order  zero-angle  iiroblem  for  f^^,  f^^^,  ...  f^^^  and  Soi  '  ’  '  ®0n  ' 

su!!H'  dolor, T.ination  must  he  made  of  Kq|,^I  •  simplest  thing  to  do 

ir,  ;  L>  truncate  t!\e  [;oj-ies  by  sotting  Ldenticaliy  equai  to  zoro, 

Tho  resulting  cqn.ations  may  then  be  integrated. 

It  is  found  that  the  difi'erentiai  equ.ition.s  fer  and  are 

itOii.  Wiit'Ifclii  u-IiU  1  Lf  i  i.  J !  1  S.  I’  d.fl  '.li'.'  T  U  “?!  f  1 1  < 

seri-'s  arc  'inear.  The  eqriticns  for  all  the  angle  of  attack  function:; 
are  ii'.o-ar  ai'.d  ccr.tair.  ccefficients  that  depend  on  the  zoro-fingle  solution. 
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Integration  is  terminated  and  the  body  shape  determined  by  the 
boundary  condition  (ll).  For  the  angle -of -attack  problem,  this  condition 
also  determines  8  . 

To  illustrate  the  procedure,  the  equations  and  boundary  conditions 
for  the  first  and  second  truncations  for  a  parabolic  shock  wave  at 
M  =  M  will  be  presented- 

After  inserting  series  (l2a)  and  (l2b)  into  equations  (8a),  (8b), 
(9b),  and  (lOb),  and  equating  coefficients  of  like  powers  of  i  and  e  , 
the  following  differential  equations  and  initial  conditions  are  obtained- 
(Primes  indicate  differentiation  with  respect  to  q  ,  The  powers  of  | 
and  €  fi'ora  whose  coefficients  the  equation  results  are  indicated  in 
front  of  each  equation) . 
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Initial  Conditions 


a)  Zero  angle 

i° }  ••  Sqo(i)  =  1  (12=) 

I  ,  ^ 

,  e°  :  SqI*'^^  "  ^ 

1^,  €°  ;  fQ^(l)  =  0  ;  f^^(l)  =  0  (l6c) 

t)  Angle  of  attack 


1°,  e  :  f^^d)  =  i  +  6  ;  f ‘^d)  =  ^  (l3d) 

I  ,  €  :  g^Q(l)  =  0  (li;d) 

5^,  e  ;  f^^d)  =  -  I  ;  f’^(l)  =  0  (I5d) 

e  :  g^^(l)  =  0  d6d) 


The  body  condition,  eq.(ll)^  becomes 

+  O(g^)  =  0  . 

(17a) 

However,  is  itself  a  function  of  |  and  depends  parametrically  on 

e  .  Expanding  e)  about  its  zero  angle  value  yields 
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Tl^d;  e)  =  (0  +  (5)  +  O(e^) 


Using  eq.(l8),  the  functions  of  in  eq.(l7a)  may  he  expanded  as 

follows 
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Eq.(l7a)  then  becomes 
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Further  expansion  of  q^*^^(|)  and  q^^^(l)  about  |  =  0  yields 


,(o) 


■b 


0  2: 


I  +  q^^^"  +  0(1^ 


0  2: 


(19) 


(20) 


-  18  - 


where 


-  19  - 


(0)  '2  (0) 


+  "o"’  “'il 


_(0)'  (0)"j,„  /„(0)\  1„(0)'3,,„(',(0) 

^  ^0  ^0  ^10^0  y  3 


^0  10 


Since  this  equation  must  hold  for  all  |  and  e  ,  equating  coefficients 
of  like  powers  of  |  and  e  yields 
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(21a) 
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Eq. (l4e)  is  the  condition  used  to  terminate  the  integration  in  any 

I'O) 

truncation  and  determines  ,  the  value  of  q  at  the  nose  of  the 

zero-angle  body.  Since  is  thus  determined,  eq.(l3e)  is  the 

condition  that  determines  5  ,  i.e. ,  this  condition  is  satisfied  for 
only  one  value  of  5  in  the  initial-condition  eq.(l3d)  and  differential 
eqs.(l3b),  (l4b),  (l5b),  and  (l6b). 


Eqs.(21a)  -  (21e),  along  with  eqs. (l8) ,  (19),  and  (20),  determine 


the  body  shape.  Eq. (21a)  yields  q 


(0) 


This  indicates  that  the 


zero -angle  body  is  represented  by  an  even  series  in  powers  of  |  .  The 

same  conclusion  could  have  been  reached  from  symmetry  cons Iderai ions. 

/  \ 

Eq.  (21b)  yields  -  O  ,  i  .e,  ,  there  is  no  perturbation  to  the  zero- 

angle  body  at  its  nose  to  this  order  of  analysis.  Eqs.(21c)  and  (21d) 
(0)"  (1)' 

determine  >  ’’to  »  I*’  terms  of  the  values  of  the  functions  of 

U  U  /■  /V  \ 

q  and  their  derivatives  evaluated  at  q^  '  .  Eq. (21e)  yields 

(i)”  ^ 

qQ  =  0  ,  which  indicates  that  the  angle -of -attack  perturbation  to 

the  body  is  represented  by  an  odd  series  in  powers  of  |  . 


Thus,  it  is  seen  that  the  condition  '^^(’^•]^)  =  0  determines  not 
only  the  body  shape,  but  also  the  unknown  quantity  S  . 

Eqs.(13a)  -  (l3e) ,  and  (l4a)  -  (l4e)  comprise  the  first-truncation 
problem,  Eqs,(l3a),  (l4a),  (l3c),  (l4c),  and  (l4e)  define  the  zero- 
angle  problem,  whereas  eqs,(l3b),  (l4b),  (l3d),  (l4d),  and  (l3e)  define 
the  angle  of  attack  problem.  Note  that  the  zero-angle  differential 
equations  (l3a)  and  (l4a)  are  nonlinear  and  that  eq. (l4a)  contains  , 

as  previously  mentioned.  Setting  g^^  equal  to  zero  allows  the  zero- 
angle  jiroblem  to  be  solved.  Setting  f^^^,  f^^  ,  and  g^^^  equal  to 
zero  in  eqs,(l3b)  and  (l4b),  and  utilizing  the  zero-angle  results, 
initial  condition  eq6,(13d)  and  (l4d),  and  body -condition  cq.(l3e) 
allows  solution  oi'  the  first-truncation  equations  for  the  angle-of- 
attack  problem. 


Similar  remarks  apply  to  the  second-truncation  problem  defined  by 
eqs.(l3a)  -  (l3e),  (l4a)  -  (l4e),  (l5a)  -  (l5e),  and  (l6a)  -  (l6e),  with 
functions  having  a  second  subscript  of  1  being  retained  and  those 
having  a  second  subscript  of  2  being  set  equal  to  zero.  Higher-order 
problems  are  solved  in  like  manner. 


2.  Three-Dimensional  Flow 

In  contrast  to  the  case  of  plane  flow,  the  continuity  equation  is 
satisfied  hy  ‘he  introduction  of  a  pair  of  stream  functions,  and 

according  tc 

d 
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pq  =  X  V\t,'2 


(22) 


The  use  of  such  a  pair  of  functions  in  this  manner  was  introduced  by 
Clebsch  [9]- 

The  entropy -conservation  equation,  (5e),  then  states  that 

^  (23) 

P 

where  the  functional  form  of  g  is  to  be  determined  from  the  boundary 
conditions  at  the  shock  wave. 

Using  equations  (22)  and  (23)  in  equations  (51>)^  (5c),  and  (5d) 
yields  tiie  following  equations  for  P  • 
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^,2  ,2  (>+l) 

+  p  ‘ 
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=  0 
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Before  the  initial  conditions  are  presented,  a  few  words  are  in 
order  regardln;;  the  stream  functions  \ifj^  and  .  These  functions 
possess  a  degree  of  arbitrariness  in  that  it  is  only  required  that  the 
surfaces  =  constant  and  =  constant  be  stream 

surfaces,  that  is,  surfaces  in  which  streamlines  are  Imbedded.  The 
intersections  of  these  surfaces  are  then  streamlines.  If  one  considers 
a  uniform  stream  for  simplicity,  there  are  any  number  of  pairs  of 
surfaces  which  can  be  set  up  to  satisfy  this  requirement.  For  example, 
one  such  pair  consists  of  orthogonal  or  shewed  planes.  Another  is 
concentric  circular  cylinders  and  orthogonal  planes  through  the  common 
axis  of  the  cylinders.  If,  however,  one  considers  axisymraetric  flow  in 
a  coordinate  system  and  chooses  for  \|/^  the  azimuthal  an;ular 

coordinate  ip  ,  \|r^  is  the  familiar  Stokes  stream  function.  Hence, 

since  we  are  considering  small  angle  of  attack,  it  is  reasonable  to 
perturb  bg  about  <p  and  about  the  Stokes  function. 

Consider  Fig. 2.  Depicted  there  is  a  shock  wave  that  is  a  paraboloid 
of  revolution  at  angle  e  to  the  uniform  free  stream.  A  line  at  aiif'le 
c  through  a  known  point  on  the  x-axis  (e.g.  ,  central  point  of  a  hyjjer- 
boloid,  focus  of  a  paraboloid,  center  of  a  sphere,  point  of  intersection 
L'f  major  and  minor  axes  of  an  ellipsoid)  is  chosen  as  reference  from 
wl'.ich  the  distance  to  the  body  streamline  ahead  of  the  shock  is  measured. 
Concentric  circular  cylinders  about  =  0  and  planes  of  angle 
P  =  constant  ai'e  cliusen  for  the  stream  surfaces  ahead  of  the  shock 
(see  section  A-/\,  Fig. 2).  V/ith  \i(,  and  defined  in  this  manner, 

the  initial  conditions  are 
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For  the  special  case  of  a  paraboloidal  shock  (B  =  0,  C  -  1)  the  stiil't, 
in  reference  line  previously  mentioned  on  pa/'e  9  results  in 
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The  boundary  condition  at  the  body  is 


i  =  0  at  T|  =  r). 


(27) 


The  method  of  solution  is  essentially  the  same  as  for  plane  flow. 
Separate  variables  by  setting 


2^1  =  i’oo(7)?^  +  f0,(r,)e  +  4-  0(5  ) 
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^  p  =  goo(ii)  +  eQ2(^)^^+ 


+  e  cos  <p 


=10' 


(n)i  +  g;^;,^(n)^  +  g-io(n)^+  o(?') 


12' 


+  oCc*^)  (28c) 


Again,  after  one  substitutes  tb.ese  expansions  into  eqs.(24)  -  (2?)  and 
equates  coefficients  of  like  powers  of  |  and  e  ,  the  series  must  be 
truncated  to  obtain  a  solution.  As  in  the  plane-flow  case,  the  equations 
for  f^Q  and  are  nonlinear  and  succeeding  equations  for  the  other 

unknown  functions  of  q  are  linear,  with  the  equations  for  the  aiigle-of- 
attack  functions  having  variable  coefficients  depending  on  the  zero- 
angle  Eolution.  The  body  boundary  condition,  eq.(27),  is  applied  in  a 
similar  manner  as  outlined  for  plane  flow  with  the  modification  that  th.c 
body  shape  is  now  given  by 

n^(t,4';  (0  e  cos</)n^^^  (I) 


3.  flonlc  Lines  and  Surface  Pressure  Distributions 

In  addition  to  body  shapes,  sonic  lines  and  surface  presstire  dis¬ 
tributions  out  to  the  sonic  line  were  calculated  to  Investigate  conver¬ 
gence  of  the  series  expansions  away  from  the  axis.  The  mathematical 
details  of  obtaining  these  quantities  will  now  be  presented. 

The  energy -conservation  equation  for  a  perfect  gas, 


y  P  4-  =  -L-  ^  4  ^ 

7-1  o  2  “  7-1  2 


(29) 


where  barred  quantities  represent  dl.mens  ionai  physical 


q uant.ities  , 


becomes,  after  suitable  nondlraensionallzation 


m2  =  2-  fg  -*•  (y  -  1)M^  P  .  ■ 

“  >■  -  H  2rM2  P  J 


(30) 


where  is  the  local  Mach  number  of  the  flow  in  the  shock  layer.  On 

the  sonic  line,  the  local  I'fach  number  is  unity.  Furthermore,  the  mass- 
and  entropy-conservation  equations  yield  p  =  f('l')  (plane  flow),  and 

P  =  o^  (three-dimensional  flow).  Putting  this  Information 

into  eq. (30)  yields 
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2  +  (y  -  l)^^  ^ 


y  +  1 
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(31b) 


Introduction  of  f  or  g  in  terms  of  the  series  expansions  for  \1/,  i|(-j , 
and  yields  equations  which,  for  a  given  truncation,  are  algebraic 

expressions  for  |  in  terms  of  functions  of  n  .  Oelection  of  values 
of  1)  between  the  shock  and  body  allows  corresponding  values  of  |  to 
be  calculated,  thus  determining  points  on  the  sonic  lino. 

The  equations  for  surface  pressure  distributions  are  obtained  by 
utilizing 


p  =  p’'  f(t)  (v  =  0)  (7) 

r  =  p’'  ^^,^2)  (v  -  1)  .  (23) 

Again,  substituting  appropriate  expansions  I'or  p,  f  ,  and  g  yields 
r(^,T;  e)  ■  t'se  of  known  pairs  of  values  of  ^  and  on  the  body 

allows  the  pressure  to  be  calculated  in  any  given  truncation. 
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III.  RESULTS  AND  DISCUSSION 


A.  Zero  Angle  of  Attack 

To  examine  the  convergence  of  the  series  expansions  for  the  stream 
functions  and  density,  zero -angle -of -attack  resttits  will  be  presented 
initially.  From  consideration  of  the  simplifications  brought  about 
in  the  equations  and  boundary  conditions,  the  free  stream  Mach  number 
is  taken  to  be  infinite  and  the  shock  shape  to  be  a  portion  of  a  circle, 
parabola,  sphere  and  paraboloid  of  revolution.  The  specific  heat  ratio, 
y  ,  is  taken  to  be  l.U. 

The  first  four  truncations  have  been  carried  out  for  the  circular 
and  paraboloidal  shock^  while  the  first  three  truncations  have  been  cal¬ 
culated  for  the  spherical  and  parabolic  shock  waves.  The  numerical 
integration  of  the  ordinary  differential  equations  for  the  unknown 
functions  of  r|  was  accomplished  using  a  fourth-order  Runge -Kutta-Gill 
method  llOJ.  The  computations  were  carried  out  on  an  IBM  7090  electronic 
computing  machine  at  Lockheed  Missiles  and  Space  Company,  Sunnyvale, 
California.  Machine  times  of  two  to  three  minutes  were  required  for 
each  case. 

Results  on  and  away  from  the  axis  are  compared  with  numerical 
solutions  obtained  using  the  inverse  method  developed  by  Van  Dyke  [l] 
and  refined  by  Fuller  [ll].  Convergence  on  the  axis  is  extremely  rapid, 
three  truncations  giving  the  stagnation -point  density  correct  to  five 
significant  figures  and  the  standoff  distance  to  four  significant  figures 
as  determined  by  comparison  of  the  parabolic  shock  solution  with 
Garabedian's  results  [12].  Convergence  off  the  axis  is  somewhat  less 
rapid,  as  is  to  be  expected  from  the  form  of  the  expansions.  The  term 
convergence  is  used  here  in  the  sense  that,  whenever  two  successive 
truncations  yield  values  for  a  physical  quantity  identical  in  three  or 
more  significant  figiures,  it  is  said  that  the  quantity  has  converged. 
Nothing  can  be  said  about  convergence  of  the  assumed  series  in  the  strict 
mathematical  sense,  since  the  functional  form  of  the  series  coefficients 
is  not  explicitly  known. 
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1.  Plane  Flow 


Results  for  the  circular  and  parabolic  shock  waves  are  given  in 
Figures  3-6  and  11.  Figures  3-5  apply  to  the  circular  shock.  In 
Figure  3^  the  convergence  of  the  reduced  stream  function  >l/s  ,  and 
the  density  function,  (>-l)p/ (y+l)  ,  along  the  axis  (  |  =  0  )  is 
examined  and  results  of  one,  two,  three,  and  four  truncations  compared 
with  numerical  solutions.  As  previously  mentioned,  convergence  is 
extremely  rapid  for  both  the  density  and  stream  function,  the  differ¬ 
ences  between  the  second,  third  and  fourth  truncations  being  indiscern¬ 
ible  to  the  scale  of  the  plot.  The  standoff  distance  referred  to  the 
shock  nose  radius,  determined  to  be  O.1895  by  the  second 

truncation  and  O.I899  by  the  third  and  fourth  truncations.  This  con¬ 
verged  result  is  smaller  than  the  value  of  0.1912  determined  by  NASA 
case  MDK* ,  and  larger  than  the  value  of  0.l88u  determined  by  NASA  case 
3iv  (these  values  for  standoff  distance  are  obtained  using  that  value 
of  ri  at  which  f  goes  to  zero,  along  with  eq,.  (2d)).  The  method 
of  solution  for  NASA  case  MDF  contains  a  smoothing  process  [11]  whereas 
that  for  case  Zk  does  not.  As  pointed  out  in  ref.  [ll] ,  tlie  smoothing 
process  has  the  effect  of  increasing  standoff  distance.  Comparison  of 
the  two  NASA  results  with  the  converged  value  leads  to  the  conclusion 
that  the  NASA  data  are  not  sufficiently  accurate  to  be  used  as  an 
absolute  basis  of  comparison.  This  point  will  be  further  borne  out  as 
other  physical  quantities  are  considered. 

Figure  U  compares  sonic  lines  and  body  shapes  as  determined  from 
the  second,  third,  and  fourth  truncations  with  the  results  ol'  NASA  case 
MDF.  The  sonic  line  determined  by  the  first  truncation  is  a  circle 
having  a  radius  of  1.593  and  center  at  x  =  1,  r  =  0  ,  while  the  first, 
truncation  body  determined  by  the  numerical  integration  is  the  point 
X  =  L,  r  =  0  .  Note  that  the  first-truncation  sonic  line  lies  ahead 
of  the  shock.  This  result  is  obtained  by  analytic  continuation  of  the 
solution  upstream  of  the  shock  wave  and  is,  of  course,  physically 

*  NA.hA  results  for  cases  flDF,  MDD,  RSI  and  RSU  courtesy  of  Ames 
Aeronautical.  I,aboratory  (\inpiubll shed) . 
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unrealizable  since  the  flow  ahead  of  the  shock  is  uniform  at  infinite 
Mach  number.  Comparison  of  the  body  shapes  determined  by  the  third  and 
fourth  truncations  indicates  convergence  beyond  the  sonic  point  to  a 
shape  other  than  the  NASA  result.  The  sonic  line,  however,  has  not 
converged  and  higher  truncations  would  be  required  to  determine  its 
exact  location. 

Figure  5  compares  the  surface  pressure  distribution  determined  by 
the  second,  third,  and  fourth  truncations  with  the  NASA  results.  Con¬ 
vergence  has  been  obtained  for  about  two-thirds  of  the  distance  to  the 
sonic  point.  The  converged  portion  of  the  curve  and  the  NASA  results 
are  parallel,  but  do  not  coincide  due  to  the  discrepancy  in  body  position. 
Higher  truncations  are  required  to  obtain  convergence  all  the  way  to  the 
sonic  point,  but  the  trend  indicates  that  the  converged  values  will 
probably  not  parallel  the  NASA  data  for  the  entire  distance. 

Figures  6  and  11  apply  to  a  parabolic  shock.  The  result.-?  of  one, 
two,  and  three  truncations  for  the  reduced  stream  function  and  density 
function  along  the  axis  arc  compared  with  both  NASA  results  and 
Garabcdlan's  solution  [12]  In  Figure  6.  Results  for  the  circular  shock 
have  indicated  that  the  NASA  data  lack  sufficient  accuracy  to  be  used 
as  an  absolute  basic  for  comparison  and  that  the  method  under  consider¬ 
ation  poGsesGCS  internal  consistency.  The  question  of  convergence  to 
tl'.c  correct  solution  may  be  answered  by  comparison  with  Garabedian's 
result  for  a  parabolic  shock  at  infinite  free  strea.m  Mach  number. 
Garabedlan  claims  l/lO  of  one  percent  accuracy  for  his  standoff  distance 
of  0. iBtlG.  The  third  truncation  yields  a  standoff  distance  of  O.1851U, 
thu:?  indicating  at  least  three  and  perhaps  four  significant  figure 
accuracy  I'or  this  quantity. 

Figure  11  compares  body  shape  and  sonic  lines  as  determined  by 
one,  two,  and  three  truncations  with  NASA  results.  For  this  case,  com¬ 
parison  indicates  near  ccnvcrge.nce  of  the  third  truncation  for  both 
quantities . 
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2.  Axisyniinetrlc  Flow 


Resultc  for  spherical  and  paraboloidal  shock  waves  are  given  in 

Figures  7  -  10«  Figure  7  depicts  the  variation  of  the  reduced  stream 
.  2 

function^  2i|//|  ,  and  density  function  along  the  axis  of  symmetry 

between  a  spherical  shock  and  the  body.  Again,  the  convergence  is 
extremely  rapid,  three  truncations  giving  a  standoff  distance,  > 

of  0.0998  compared  with  the  NASA  results  of  O.O98I  (unsmoothed)  and 
0.0993  (smoothed). 

Figures  8-10  apply  to  a  paraboloidal  shock  wave.  Figure  8  yields 
a  standoff  distance  of  O.8956  as  determined  by  the  second,  third,  and 
fourth  truncations..  Hence,  the  standoff  distance  converges  to  four 
significant  figures  in  only  two  truncations  for  this  case. 

Figure  9  compares  body  shaxjes  and  sonic  lines  with  a  MSA  numerical, 
solution.  In  the  first  truncation,  both  body  and  sonic  lines  are 
r;  =  constant  lines,  the  body  being  at  t|  =  0.8552  and  the  sonic  line 
at  r)  -  1,0780.  The  corresponding  values  of  x  and  r  may  be  obtn.incd 
from  eqs,(2b)  and  (2c).  Again,  the  first-truncation  sonic  line  lies 
ahead  of  the  shock,  this  result  having  been  obtained  by  analytic  con¬ 
tinuation  of  the  solution  upstream  of  the  shock  wave.  Since  the  sonic 
point  on  the  body  has  a  5  coordinate  greater  than  in  the  case  of  the 
circular  shock,  neither  the  body  shape  nor  the  sonic  line  ap])ears  con¬ 
verged  even  in  four  truncations  for  this  case.  However,  consider  the 
surface  pressure  distribution  given  in  Figure  10.  One -diinensional 
streamtube  relations  yield  the  result  that,  for  7  =  1 . U  ,  the  pressure 
at  the  sonic  I'olnt  on  the  body  must  have  the  value  0.48<j.  This  vaiue  i 
plotted  at  the  value  of  x  determltsed  by  the  intersection  of  the  fourth- 
truncation  body  and  sonic  line  from  Figure  9  (asterisk  symbol ,  Figure 
10).  One  sees  that  this  point  nearly  coincides  with  the  pressure  curve 
obtained  by  applying  a  Shank's  transformation  [13]  to  the  fourth- 
truiicnticn  pressure  series.  The  nonliiiear  Shank's  transforrnatitjn 
causes  the  slowly  convergent  pressure  series  to  converge  more  rapidly. 

The  above  agreement  Indicates  near  convergence  without  the  nece.s.sity 
of  c.'ilculatir.g  the  fifth  truncation. 
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B.  Small  Angle  of  Attack 

The  analysis  outlined  for  the  angle -of -attack  problem  has  been 
carried  out  for  shock  waves  that  are  parabolas  (plane  flow)  and  parab¬ 
oloids  of  revolution  (three-dimensional  flow).  The  first  three  trun¬ 
cations  have  been  calculated  for  each  shock  shape.  The  free  stream  Mach 
number  is  taken  to  be  infinite,  and  y  takes  on  the  vclues  1.^,  .1.2, 

1.1,  1.05.  The  resulting  body  shapes,  sonic  lines,  and  positions  of  the 
body  streamlines  at  their  points  of  intersection  with  the  shock  wave  are 

shown  in  Figures  12  and  13  for  e  =  10°,  y  =  1.4,  where  a  finite  value 

2 

of  e  has  been  selected  for  graphical  purposes  and  the  error  is  0(e  )  . 
For  purpose  of  comparison,  body  shapes  and  sonic  lines  for  n  parabolic 
and  paraboloidal  shock  at  zero  angle  are  given  in  Figures  11  and  9 
respectively.  Computing  machine  times  of  five  to  six  minutes  per  case 
were  required  since  the  zero-angle  functions  appear  as  coefficients  in 
the  angle -of -attack  equations  and  hence  must  be  stored  in  the  machine. 
Also,  in  certain  cases  it  was  necessary  to  use  very  small  step  sizes 
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(e.g. ,  Aq  =  10  )  in  order  to  obtain  accurate  values  for  the  quantity 

6  . 

Of  particular  Interest  in  the  angle -of -attack  problem  is  whether 
or  not  the  streamline  that  wets  the  body  crosses  the  shock  normal.ly, 
thus  having  maximum  entropy.  For  this  to  be  so,  the  quantity  b  (see 
Figure  2)  must  have  the  value  g  .  The  results  of  the  analyses  are 
given  in  Table  I. 


TABLE  I 
Values  of  5 

Parabolic  Shock,  M  =  qq  Paraboloidal  Shock,  M  =  cq 

y  y 


Truncation 

1,4 

1.2 

1.  ] 

1.05 

l,h 

1.2 

1.1 

1.05 

1518 

8417 

- 

- 

88.00 

37.48 

25.79 

-1.34 

2 

.1.833 

.4933 

.4S7,'4 

.4990 

.  5576 

.5394 

.51 '90 

.  5042 

"7 

.47^3 

.4903 

.49(;6 

.4988 

.4871 

.4950 

.4984 

.  500it 

For  all  cases  the  trend  of  decreasing  6  with  increasing  trunca¬ 
tion  is  toward  a  value  of  ?i  slightly  less  than  4  .  For  the  three- 
dimensional  case,  a  reversal  of  this  trend  in  higher  truncations  could 
conceivably  result  in  a  value  of  ^  .  However,  even  if  this  reversal 
took  place  in  two  dimensions,  a  con'/erged  '/alue  would  be  expected  to 
lie  between  the  second-and  third-truncation  values  and  thus  could  never 
be  ^  .  Hence,  the  streamline  that  wets  the  body  passes  through  the 
shock  wave  slightly  above  the  point  where  it  is  a  normal  shock.  However, 
as  the  Newtonian  limit  process  (  M  -»  <»,  r  -*  1  )  is  carried  out,  the 
trend  is  toward  a  value  of  ^  .  Hence  in  this  limit,  when  the  parabolic 
and  paraboloidal  shock  and  body  coincide,  the  body  is  covered  by  the 
maximum-entropy  streamline  in  both  the  plane  and  three-dimensional  cases. 

Since  the  shock  waves  possess  symmetry  about  the  x--axls,  the 
resulting  body  shapes  are  asymmetric  to  this  axis.  Consideration  of 
the  comparison  In  Figure  11  of  the  third-truncation  body  for  the  para¬ 
bolic  shock  at  zero  angle  with  the  NASA  results,  and  of  the  results  for 
the  circular  shock  given  in  Figure  4,  lead  to  the  conclusion  that,  for 
ayminel.ric  plane  flow  at  M  -  ■»  ,  three  truncations  yield  a  nearly  con¬ 
verged  body  out  to  the  sonic  point.  Since  e  is  taken  as  small,  it 
is  reasonable  to  extend  this  conclusion  to  include  small  angles  of 
incidence.  With  the  above  facts  in  mind,  and  the  realization  that  It 
is  somewhat  impractical  to  t.alk  about  asymmetric  shapes  from  an  experi¬ 
mental  or  engineering  viewpoint,  an  attempt  was  made  to  find  an  axis 
system  in  wlilch  that  portion  of  the  two-dimensional  body  given  by  the 
third-truncation  (Figure  12)  out  to  the  sonic  points  could  be  closeiy 
approximated  by  a  conic  section.  An  ellipse  was  fitted  by  a  trlal-and- 
error  procedure  in  wliich  the  positions  of  the  foci  and  ms,jor  axis  were 
varied  until  points  on  the  ellipse  and  third  truncation  body  coincided 
to  three  significant  figures  (Figure  lA).  The  ellipse  is  given  by  the 
equation  x‘"/0.37  +  r  /0.2”  =  1  ,  where  the  origin  of  the  (  x,  r  )  axes 
is  at,  X  =  .”90,  r  =  -  .0.53  ,  and  the  x  axis  is  at  angle  Ib.PO*^ 
measured  clockwise  from  the  free  stream  direction.  Thus,  the  shock 
wave  is  seen  to  rotate  about  '0  per  cent  as  fast  as  the  body  as  t'te 
angle  of  attack  is  Imposed. 


-  33 


No  attempt  was  made  to  perform  a  corresponding  fit  to  the  three- 
dimensional  bod,y  determined  by  the  third  truncation,  since  Figure  9 
indicates  that  the  body  shape  out  to  the  sonic  point  has  not  converged. 
Note,  however,  that  the  third-truncation  body  is  closer  to  the  shock  in 
the  cp  =  TT  plane  than  in  the  (p  -  0  plane.  Hence,  to  this  truncation, 
the  shock  wave  still  .must  rotate  more  slowly  than  the  body  as  angle  of 
attack  is  Imposed. 
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IV.  CONCLUSIONS 


A  semi -analyt leal j  semi -numerical  method  of  analysis  has  heen 
developed  and  successfully  applied  to  the  problem  of  two-and  three- 
dimensional  hypersonic  flow  past  blunt  bodies  at  small  angles  of 
attack. 

Solutions  at  zero  angle  of  attack  are  obtained  as  a  special  case 
of  the  general  problem.  Standoff  distances  are  predicted  to  higher 
accuracy  than  by  any  other  existing  method. 

For  a  shock  wave  that  is  a  parabola  or  paraboloid  of  revolution 
at  angle  of  incidence  to  a  uniform  free  stream  having  infinite  Mach 
number  and  y  =  1.4  ,  the  streamline  that  wets  the  body  passes  through 
the  shock  wave  slightly  above  the  point  wliere  it  is  a  normal  shock  wave. 
Hence,  the  body  is  not  covered  by  the  maximum-entropy  streanil.lne  in 
these  eases.  In  the  Newtonian  limit  (  M  =  «>  ^  y  =  1  ) ,  the  parabolic 
or  paraboloidal  body  is  wetted  by  the  maximum-entropy  streamline.  For 
a  parabolic  shock  at  infinite  Mach  number,  7  =  1.4  ,  the  shock  wave 
rotates  more  slowly  than  the  body  as  angle  of  attack  is  Imposed. 
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V.  FUTURE  RESEARCH 


Although  T,he  analysis  presented  here  is  for  general  Mach  number 
and  shock  shape,  results  have  been  obtained  only  for  certain  shock 
shapes  at  infinite  Mach  number.  Work  is  currently  in  progress  to 
program  the  general  problem  for  the  IBM  7090  computer. 

Since  the  results  of  these  calculations  in  the  vicinity  of  the 
sonic  line  might  serve  as  input  data  to  a  characteristic  program  for 
the  solution  of  the  flow  field  downstream  of  the  sonic  line,  possible 
means  of  improving  the. convergence  of  the  method  in  the  vicinity  of  the 
sonic  line  are  being  explored. 

One  such  possibility  involves  expanding  the  stream  functions  and 
density  in  Taylor  series  about  an  arbitrary  |  =  constant  line  in  the 
flow.  By  choosing:  this  arbitrary  line  close  to  the  sonic  line,  one 
might  expect  more  rapid  convergence  for  values  in  this  region  at.  a 
sacrifice  of  convergence  in  the  vicinity  of  the  shock  axis  of  symmetry. 

Another  possible  means  of  improving  convergence  away  from  the 
shock  axis  of  symmetry  is  by  applying  power  series  expansions  to  quan- 
tit.ies  wlilch  vary  slowly  as  one  moves  away  from  the  shock  symmetry  axis, 
thereby  reducing  the  magnitude  of  the  coefficients  of  the  higher  order 
terms  with  respect  to  those  of  lower  order.  For  example,  at  infinite 
!4nch  number,  p  is  constant  behind  the  shock,  whereas  p/p^  is 
constant  on  the  body.  If  one  considers  expandinf:  the  product  p™  p^  , 
then  m  -  1  corresponds  to  expanding,  p  ,  whereas  m  -  T'/t'-I  corresponds 
to  expanding  p/m^  .  One  migiit  consider  expanding  for-  a  value  of  m 
wliich  is  the  average  of  1  and  y/y-\  .  This  would  then  correspond  i,o 
a  quant,  ily  which  varies  more  slowly  as  one  moves  away  frotn  the  shock 
axis  of  symmetry  than  either  p/p^  on  the  s'-'ock  or  p  on  the  body. 

A  further  possible  extension  cf  the  work  is  to  render  it  applicable 
to  hypersonic  flow  of  real  gases  in  equilibrium  and  perfect  yj'ases  in 
t.beriKi ,1  and  ch.entical  nonequ llibr ium. 
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Fiji.  1  Representation  of  Shock  Wave 
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Fig.  2  Parabolic  or  Paraboloidal  Shock  Wave  at  Angle  of  Attack 


Fig.  3  Variation  of  Reduced  Stream  Function  <li/i  and  Density  Along  the  Axis 

Circular  Shock,  M  =  "  ,  y  =  1-  4 


Fig.  5  Surface  Pressure  Distribution  —  Circular  Shock, 
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Fig.  6  Variation  of  Reduced  Stream  Function  ij/i,  and  Density  Along  the  Axis 

Parabolic  Shock,  M  =  *  ,  7  =  1-4 


Fig.  7  Variation  of  Reduced  Stream  Function  2i/£  and  Density  Along  the  Aicis 

Spherical  Shock,  M  =  *  ,  7  =  1-4 
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Variation  of  Reduced  Stream  Function  24-/^  and  Density  Along  the  Axis 
Paraboloidal  Shock,  M  =  ,  y  =  1.4 


Fi^;;.  !)  Comparison  of  Body  Shapes  and  Sonic  Lines 
Paraiioloidal  Shock  ,  M  ^  =  1 . 4 


■iH 
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Fig.  10  Surface  Pressure  Distribution  —  Paraboloidal  Shock,  M 


Fig.  13  Body  Shapes,  Sonic  Lines,  and  Comparison  of  Body  Streamlines  With  Maxi¬ 
mum  Entropy  Streamline  —Paraboloidal  Shock,  M  =  'w,y  =  1.4,e  =  10° 
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Fitting  an  Ellipse  to  tlie  Third  Truncati?  ?|  v 

M  =  "o  .  y  =  1,  //  .- 
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